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, Dimensional Reduction is applied to QCD in order to compute various renor- 

malization constants in the DR scheme at higher orders in perturbation theory. In 
particular, the j3 function and the anomalous dimension of the quark masses are 
derived to three- loop order. Special emphasis is put on the proper treatment of the 
so-called e-scalars and the additional couplings which have to be considered. 
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& 1 : 1 Introduction 

si ■ 

The big success of the Standard Model of elementary particle physics is based to a large 
extent on precision calculations which sometimes reach three-, four- and even five- loop 
accuracy. Such calculations currently all rely on Dimensional Regularization (DREG) [1, 
2] which is an elegant and powerful tool to parametrize the divergences occurring at 
intermediate steps of the calculations. 

In DREG, the number of space-time dimensions is altered from four to D = 4 — 2e, 
which renders the loop integrations finite. It is clear, however, that if DREG is applied 
to a 4-dimensional supersymmetric theory, the number of bosonic and fermionic degrees 
of freedom in super-multiplets is no longer equal, such that supersymmetry (SUSY) is 
explicitly broken. In order to avoid this problem, Dimensional Reduction (DRED) has 
been suggested as an alternative regularization method [3]. Space-time is compactified to 
D = 4 — 2e dimensions in DRED, such that the number of vector field components remains 
equal to four. Momentum integrations are D-dimensional, however, and divergences are 
parametrized in terms of 1/e poles, just like in DREG. Since it is assumed that e > 0, 
the four-dimensional vector fields can be decomposed in terms of D-dimensional ones 
plus so-called e-scalars. The occurrence of these e-scalars is therefore the only difference 
between DREG and DRED, so that all the calculational techniques developed for DREG 
are applicable also in DRED. 



Nevertheless, it soon was realized that DRED suffers from mathematical inconsistencies 
in its original formulation [4]. Currently, it seems that they can only be avoided by 
interpreting the fields as living in an infinite dimensional space, which again leads to 
explicit SUSY breaking [5,6]. A higher order calculation will therefore require similar 
SUSY restoring counter terms as they are needed in DREG in general. For some of the 
currently available two-loop results, however, it has been shown that these counter terms 
vanish [7]. 

Although DRED was originally constructed for applications in supersymmetric mod- 
els, it has been shown that in certain cases it can be useful also in non-supersymmetric 
theories [8-10] like QCD. For example, since it is possible to turn QCD (with massless 
quarks) into a super- Yang-Mills theory by simply adjusting the colour factors, a calcula- 
tion using DRED provides the possibility to use non-trivial Ward identities for a check of 
complicated calculations (see, e.g., Ref. [11]). 

As mentioned before, DRED parametrizes ultra-violet divergences as poles in e. One 
can therefore formulate a renormalization scheme analogous to the MS scheme, usually 
called the DR scheme. In this paper, we compute the beta function of the strong coupling 
and the anomalous dimension of the quark masses to three-loop accuracy within this 
scheme. An important issue turns out to be the renormalization of the qqe vertex. It 
requires to introduce a new, so-called evanescent coupling constant a e . A similar argument 
holds for the four-e-scalar vertex, but at the order considered here, this vertex does not 
get renormalized. The proper treatment of a e leads us to conclude that the three-loop 
result for the QCD (3 function available in the literature [12] is incorrect. The correct 
result is provided in Section El 

The outline of the paper is as follows. In Section El we provide the notation and set 
the general framework for the calculation. Subsequently, we describe in Sections EH and H] 
the calculation for the (3 and the 7 m function up to three loops. Section contains the 
conclusions. 

2 Framework 

We consider QCD and apply Dimensional Reduction (DRED) as the regularization scheme. 
Thus, besides the usual QCD Feynman rules for quarks (q) and gluons (g), we have to 
consider additional vertices involving the so-called e-scalars, namely qqe, gee, ggee, eeee 
(for the corresponding Lagrange density, see, e.g., Refs. [13,14]). In general, also a mass 
term for the e-scalar has to be taken into account [15,16]. However, on simple dimensional 
grounds it affects neither the QCD /3-function nor the anomalous dimension of the quark 
mass, so we do not need to consider it here. 

In a non-supersymmetric theory, it is important to note that the qqe and the qqg 
vertices renormalize differently. Therefore, one needs to distinguish the coupling constant 
g e , multiplying the qqe vertex, from the strong coupling g s [8]. Also the eeee vertex 
renormalizes differently; in fact, in QCD one needs to allow for a more general colour 
structure of this vertex, leading to three additional coupling constants X r (r — 1, 2, 3). In 
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order to fix the notation we display the relevant part of the Lagrange density [8] 

C = ...-^\ r H? bcd ele c a ,e b a e d a , + ... : (1) 

r=l 

where e denotes the e-scalar fields, and a and a 1 are 2e-dimensional indices. For the SU (3) 
gauge group, the Hf >cd are three independent rank four tensors which are symmetric under 
the interchange of (ab) and (cd). Our choice 

jjabcd ^_ ^ j-ace j-bde _|_ jade j'^ ce ^ 

jjabcd ^ab^cd _|_ ^ac^bd _|_ ^ad^bc 

jjabcd = ( S ac 6 bd + 5 ad 5 bcj _ fiabficd ^ ^ 

fixes the Feynman rules in a unique way. Note that for SU(N C ), N c > 3, there are four 
independent tensors Hf' cd . At the order considered in this paper no renormalization 
constant for the eeee vertex has to be introduced. The vertices gee and ggee, on the 
other hand, are renormalized according to g s because of gauge invariance [8]. 
g e and A r will be called "evanescent couplings" in what follows, and we define 

9s 9e i ""V / Q x 

<Xs = -7-, «e = — and i] r . = — . (3) 

47T 47T ATT 

The renormalization constants for the couplings g s and g e , the quark mass to, the QCD 
gauge parameter £, as well as for the fields and the vertices are introduced as 

9° s = V £ Z s g s , g° e = /i e Z e g e , m° = mZ n 

l-e°=(l-0^3, q°=Vz~2q, G^ = y/Z~ S G^, 
e°f = ^Zje a a) 

^qqG = Zl^qqG , ^qq £ = Z[I \ qe , T%q = Z X Y c5G , (4) 

where \x is the renormalization scale, D = 4 — 2e is the number of space-time dimensions, 
and the bare quantities are marked by the superscript "0". The quark, gluon, e-scalar, 
and ghost fields are denoted by q, G^, e a and c a , respectively, and T xyz stands for the 
vertex functions involving the particles x, y and z (a is the colour index.). The gauge 
parameter £ is defined through the gluon propagator, 



9l 


= I^Z e g e , 


m° 




= \fZ2q, 




c 0,a 


= \/z^c a , 




r°- 


— 7 £ V 

— ^l 1 qqe i 


1 ccG 



9 ^-i 
q 2 + ie 



From the renormalization of the ghost-gluon or quark-gluon vertex one obtains the 
renormalization constant of the strong coupling 
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Similarly, the quark-e- scalar vertex leads to the relation 



7 e 

Z e = — (7) 

Z2 v/ Z3 

It is well-known that Z s 7^ Z e even at one-loop order. Furthermore, both Z s and Z e 
depend on g s , g e , and X r [8]; note, however, that Z s depends on g e and A r only starting 
from three- and four-loop order, respectively, while Z e depends on g e and A r already at 
one- and two-loop order, respectively. 

Let us next introduce the (3 functions both for DREG and DRED. In DREG, of course, 
the e-scalars are absent, and from the definition 

P («. ) - A* ^ — (8) 

the usual relation between /3 MS and Z s is obtained: 

^ms / « 1ti 7 ms' 

(3 MS { MS) = _ e ^_ 1+2a MS^| (Q) 

where denotes Z s evaluated in the MS scheme, and is the usual definition of the 
strong coupling within this scheme [17]. /3 is known to four-loop order (see Refs. [18,19] 
and references therein). Due to the fact that we have five different couplings in DRED, 
the relations between Z s and Z e and the corresponding beta functions are slightly more 
involved. They are given by 

/f>FU,R.}) = ^ 2 d a?R 

d/T 7T 

« e „a e dZ e ^fyp « e <9Z e n \ / n a e dZ e \ 1 

e T + 2 *a^ + 2 *£a^)( 1 + 2 *^J ' (10) 

where it is understood that the renormalization constants Z e and in Eq. (|1U|) are 
evaluated within DRED with (modified) minimal subtraction, and a° R is the correspond- 
ing strong coupling constant in this scheme. As in the MS scheme, the coefficients of the 
single poles fully determine the f3 functions. Let us remark that the terms proportional 
to P Vr , the beta functions corresponding to the couplings r] r , contribute to /3° R only at 
the four- loop order. Furthermore, only the approximation /3 Vr = — is needed for the 
two- loop calculation of (3 e . 
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In analogy to Eqs. Q and (jlUJI we introduce the anomalous mass dimensions which 
are given by 



,.i j aiv, 7MS 

^ms^ms^ = _/f — m MS = tt3 ms — 

d/i 2 M cta^s 



m 



MS 



ni 



d/i 2 



= ~da^~ ~ Pe ~d^ e — ^^-ejT- 

As in the case of the (3 function, 7° R also gets additional terms due to the dependence of 
Z m on the evanescent coupling g e and on the quartic ^-scalar couplings A r . The four-loop 
result for 7^ s can be found in Refs. [20,21]. 

Let us add a few remarks concerning the meaning of the evanescent coupling a e at 
this point. In a non-supersymmetric theory, ot e can be set to an arbitrary value a e at an 
arbitrary, fixed scale fi, a e (jl) = a e . This corresponds to a choice of scheme and in turn 
determines the value of a^ R through, say, an experimental measurement. At any scale /x, 
both a° R and a e are then determined by the renormalization group equations fjlOj) . One 
particular scheme choice would be to set a e (fi) = af R (/i). Note, however, that already 
at one-loop level one will have a e (/i) 7^ a^ R (/i) for any /i^/i due to the difference in the 
renormalization group functions f3 s and j3 e . 

In a supersymmetric theory, on the other hand, one necessarily has /3f R = (3 e and 
a° R = a e at all scales. Thus, if one assumes that QCD is a low energy effective theory 
of SUSY-QCD, a e is no longer a free parameter. Rather, a e and a^ R are both related 
to the unique SUSY-QCD gauge coupling by matching relations (see, e.g., Ref. [22]) and 
renormalization group equations. 

These considerations show that the choice a e = a^ R is not compatible with the renor- 
malization group evolution of these couplings unless all SUSY particles are taken into 
account in the running. In fact, it cannot be assumed at any scale as soon as one or more 
SUSY particles are integrated out. An example where this is relevant already at one- loop 
level is the m DR <-> m MS relation as will be pointed out in connection with Eq. (j22j) below. 
An analogous discussion holds also for the evanescent couplings rj r . 



3 (3 function to three-loop order 

Within the framework of DRED outlined in the previous section we have computed Z 1; 
Z2, Z3, Z\ and Z3 to three-loop order. They are obtained from the two- and three- 
point functions according to Eq. (J3J) (see, e.g., Ref. [23] for explicit formulae). Thus, 
according to Eq. Z s is computed in two different ways and complete agreement is 
found. Furthermore, we compute Z\ and Z| to two-loop order and hence obtain Z e to 
the same approximation. 
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Since only the divergent parts enter the renormalization constants, we can set all 
particle masses to zero and choose one proper external momentum in order to avoid 
infrared problems. For the generation of the about 11,000 diagrams we use QGRAF [24] 
and process the diagrams with q2e and exp [25,26] in order to map them to MINCER [27] 
which can compute massless one-, two- and three-loop propagator-type diagrams. 

The n-loop calculation leads to counter terms for g s , g e , and the gauge parameter 
£, which are then inserted into the (n + l)-loop calculation in order to subtract the 
sub- divergences. We remark that e-scalars are treated just like physical particles in this 
procedure. 

For the [3 function, to a large extent it is possible to avoid the calculation with e-scalars 
and evaluate the Feynman diagrams by applying only slight modifications as compared 
to DREG. For that, after the projectors have been applied and the traces have been 
taken in D = 4 — 2e dimensions, one sets e = 0. The evaluation of the momentum 
integrals, however, proceeds in D dimensions, just as for DREG. During the calculation it 
is necessary to keep track of the qqg vertices since the difference between DREG and DRED 
in the results of the corresponding diagrams effectively accounts for the contributions from 
the qqe vertex. Thus, the renormalization constant Z e has to be used for this contribution. 

We refrain from listing explicit results for Z s but instead present the results obtained 
from Eq. fi"0|) . Although (3 e is only needed to one-loop order for the three- loop calculation 
of /3° R we present the two-loop expression which enters the three-loop calculation of 7° R . 
Writing 




7T z — ' J \ 7T / V 7T / \ 7T / V 7T / V 7T 

i ,j,k,l,m 



/3 e (a° R , a e , R.}) = -e— - Ptjkh 



" e w ( a ° R I f a e\ j fVi\ k fm\ l fV3 



/ , i i rein l i i \ / \ / \ / \ 

7T z — ' \ 71 I \ 7T / \ 7T / \ 7T / V7T 

i,j,k,l,m 



(12) 

frj 2 \ l /"r]3\ m 
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we find for the non-vanishing coefficients up to three respectively two loops: 



where 
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(13) 



(14) 



are the usual colour factors of QCD, and n/ is the number of active quark flavours. In 
Eq. (fTzTj) we introduced five indices for the coefficients of AP R an d Ae- However, we drop 
the last three indices whenever there is no dependence on r] r . In particular, f3 s depends 
on the T] r only starting from four-loop order. Note that those terms involving r] r are only 
valid for N c = 3, whereas the remaining ones hold for a general SU(N C ) group. 

As a first check on the results given in Eq. (JXEJ), we specialize them to the supersym- 
metric Yang-Mills theory containing one Majorana fermion in the adjoint representation, 
by setting C A = C F = 2T, rif = 1, and a° R = a e = T]i and r] 2 = r] 3 = 0. Accordingly, we 
obtain for the non-vanishing coefficients of the 01 



3 m 



in agreement with Ref. [28]. Moreover, comparing these coefficients for pure QCD with the 
literature, one finds that the two-loop result for /3° R and the one-loop result for (3 e agree 
with Ref. [8]. Actually, up to this order, the result for the first two perturbative coefficients 
of j3 s is the same in the DR and the MS scheme which is a well-known consequence of 
mass-independent renormalization schemes. However, our three-loop result for 0® R differs 
in the terms proportional to CpTnf, CaCfThj and CVT 2 n 2 from the one that can be 
found in Ref. [12]. 

In order to explain this difference, let us have a closer look at the method used in 
Ref. [12]. The function /3° R was derived from the known result for /3^ s by inserting the 
relation between a° R and a^ s . The couplings a e and af R , as well as their /3-functions 
P e and j3f R were identified throughout the calculation. But as we will show shortly, this 
identification makes it impossible to obtain consistent higher order results. Keeping the 
couplings different, on the other hand, the relation between aj* R and reads 



Ti 12 \ n 72 A n n 8 } 



(16) 



where the dots denote higher orders in a^ s , a e , and r\ r . We obtained this relation by 
noting that the value of a s in a physical renormalization scheme should not depend on 
the regularization procedure: 

a f = (^.x) 2 a * , ^ h - x = Zf/Zf - x , X G {MS, DR} 

/ 7 ph,DR 7 MS\ 2 (17) 

where 2'^ IS / DR are the charge renormalization constants using minimal subtraction in 
DREG/DRED, as defined above. For £P h > MS / DR ; on the other hand, we use DREG/DRED 
combined with a physical renormalization condition. We observe that the ratio in Eq. (|17|) 
is momentum independent, such that the calculation amounts to keeping the constant 
finite pieces in the charge renormalization constants ^P h ' MS / DR # Note that the various Z s 
in Eq. (JT7J) depend on differently renormalized a s , so that the equations have to be used 
iteratively at higher orders of perturbation theory. 
Equation (fTB]) has to be inserted into 

«a, M) = = ^^^5 + ^> + • • • ' (18) 

where the first equality is due to the definition of (3f R and the second one is a consequence 
of the chain rule, with terms arising through the i] r represented by dots. Using the three- 
loop expression for j3^ s (see Refs. [18,19] and references therein), we obtain the same 
result as in Eq. (|T3|) which not only provides a powerful check on the various steps of the 
calculation, but also confirms the equivalence of DREG and DRED at this order 
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Let us stress that even if one sets a e = a° R in the final result (cf. Eq. ([T3])). one does 
not arrive at the expression for (3f R provided in Ref. [12]. 

Indeed, a way to see that the identification of a[p R and a e at intermediate steps leads 
to inconsistent results is as follows. Whereas in the case of the (3 function the error is a 
finite, gauge parameter independent term, it leads to a much more obvious problem for 
the quark mass renormalization: Z m will contain non-local terms at three-loop order if 
g e = g s is assumed throughout the calculation, and 7 m as evaluated from Eq. (fTTj) will 
not be finite. 

4 Mass anomalous dimension to three loops 

In this section we use the framework of Section El in order to obtain the anomalous 
dimension of the quark masses within DRED as defined in Eq. (jllj) . The result will be 
derived both by a direct calculation of the relevant Feynman diagrams in DRED, as well 
as indirectly by using the result from DREG and the MS-DR relation between the strong 
coupling and quark mass. 

The evaluation of Z m to three-loop order proceeds along the same lines as for the 
renormalization constants of the previous section. However, in contrast to Z s , the coupling 
a e already appears at one-loop order. Thus the two-loop expression for Z e is required 
which can be obtained from Eq. (|13|). At one- loop order we find complete agreement with 
the result given in Ref. [8] ; the two-loop term is — to our knowledge — new. 

From the three-loop result for Z m we obtain the anomalous dimension 
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(20) 



Again the last three indices are suppressed whenever there is no dependence on r\ r . Fur- 
thermore, those terms involving r\ r are only valid for N c = 3, whereas the remaining ones 
hold for a general .St/ (iV c ) group. 

On the other hand, 7° R can be derived indirectly from the MS result obtained within 
DREG. The analogous equation to (fTHJ) is given by 



7E R (a? R ,ae,-for}) = 7 



MS 



<91nm DR 7r/l MS «9m DR 7r/3 e «9m DR 

= + _ — ==• H = 

«91nm MS m DR da™ s m DR <9a e 



+ . . . 



(21) 



which requires the two-loop relation between m DR and m MS in order to obtain 7° R to 
three loops. The two-loop relation between m DR and m MS can be computed in close 
analogy to Eq. (|17j) by keeping not only the divergent but also the finite parts in the 
calculation of the fermion propagator. Our result reads 
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where the dots denote higher orders in orf , a e , and rj r . The one- and two-loop terms of 
Eq. (|2~2"J) agree with Ref. [30] in the limit a e = a° R . Let us remark that in order to get 
Eq. (|22I), also the one- loop relation between the DR and MS version of the gauge parameter 
is a necessary ingredient unless one works in Landau gauge (£ = 1). We performed the 
calculation for general covariant gauge and use the cancellation of the gauge parameter 
in the final result welcome check. 

As a further ingredient we need 7^J S which can be found in Refs. [20,21]. Inserting 
this result and Eq. (|2~2~j) into (|2*T|) leads to Eq. (j2Hj) . Again, this is a powerful check on 
our calculation and shows the equivalence of DRED and DREG at this order. Note that 
in the indirect approach the i] r enter only through the factor f3 e in Eq. (|2ip. 

The two- loop result of 7° R can also be found in Ref. [30] and we agree for a e = a s . 
The three-loop result for 7° R is new. 

The distinction between a s and a e in Eq. (}2*2*|) is essential for phenomenological 
analyses as can be seen from the following numerical example. Assuming a super- 
symmetric theory and integrating out the SUSY particles at // = M z , we may use 
a® R (M z ) = a e (M z ) = 0.120 as input, and then evolve a e and a^ R separately to lower 
scales by using Eqs. (HUJ), (fT2|) . and (fTSjl . For fi b = 4.2 GeV, we arrive at 1 af K (// 6 ) = 0.218 
and a e (/ib) = 0.167, for example. Using m^ /[S (/i b ) = 4.2 GeV, Eq. (|2*2*Jl then leads to 
mf R (fj,b) = 4.12 GeV. If one wrongly identifies a e with af R in Eq. (J22J) . one obtains a 
value for m^ R (fib) which is roughly 30MeV smaller than that. Note that this difference 
is of the same order of magnitude than the current uncertainty on the 6-quark mass 
determination (see, e.g., Ref. [31]). 

Note that the identification a e = a s has also been made in Eq. (26) of Ref. [32] for 
\i = M z (see also Ref. [33]), which incorporates our Eq. (J22J) for rif = 5. This induces an 
inconsistency of order o? s {M z \ whose numerical effect is quite small. 

5 Conclusions 

In many cases, DRED poses an attractive alternative to DREG — not only for super- 
symmetric theories. We computed the QCD renormalization group function of the strong 
coupling constant (/3) and of the quark masses (j m ) to three-loop order in this scheme 
using two different methods. The agreement of the results obtained in both ways con- 
firms the equivalence of the DR and the MS renormalization scheme at this order, in the 
sense that they are related by an analytic redefinition of the couplings and masses [29]. 
Furthermore, we find that the three-loop /5-function found in the literature differs from 
ours. We trace this difference to the fact that the evanescent coupling of the qqe vertex 
had been identified wrongly with a s in Ref. [12]. 

Let us stress that higher order calculations within the framework of DRED should also 
be useful in the context of the Minimal Supersymmetric Standard Model where precision 
calculations will be important in order to be prepared for measurements at the CERN 

1 Only one-loop running of a s and a e is applied here. 
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Large Hadron Collider and other future high energy experiments. 
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